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SUMMARY 


This volume summarizes the results of an image space analysis and optim- 
ization and of a performance analysis of the telescope assembly in the 
presence of various alignment errors. It also includes a brief study 
on possible test arrangments for 1.2m or 1.5m diameter x-ray telescopes. 
Omitted is a sianmary of the work that was done on the test flats because 
of its tentative state and because much of it was already out of date at 
the time this report was prepared. 
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1. INITIAL TELESCOPE ANALYSIS 


Design parameters fo- all six subsystems compatible with the Input 
requirements of the ray trace program KROGRAZ were generated from data 
provided by SAO to MFSC. The most Important parameters together with 
some analytical results are summarized in Table 1,1. Fig. 1.1 illustrates 
the geometric meaning of the following parameters: 

P Q1 : Center radius of primary, 

P Q2 : Center radius of secondary, 
d : Center to center distance, 
b : Back focal distance 

Y : Grazing angle at the center of each element. 


PRIMARY 



Fig. 1.1: Two-Mirror Grazing Incidence System with Design Parameters 
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A ray trace analysis was performed to first find the best focal sur- 
face of all Individual subsystems and then determine ths optimum positions 
of Image planes such that the rms spot size on axis does not exceed 0.5 
arcsec. The corresponding distances of these image planes from the 
gausslan focal point and the maximum achievable half field angles provid- 
ing a spot size of no more than 0.5 arcsec for every individual subsystem 
are included in Table 1. An illustration of the image space geometry is 
shown in Fig. 1.2. For these analyses a reflectivity of 1 was assumed 
for all surfaces. 



Fig. 1.2: Meridional Section of Image Surface 

(R = Radius of Image Curvature) 
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TABLE 1.1: System Parameters 


SUBSYSTEM 


1 

2 

3 

4 

5 

6 

Center Radius of First Mirror 

(in. ) 

23.88 

21.53 

19.23 

16.93 

14.78 

12.63 

Center Radius of Second Mirror 

(in.) 

22.8273 

20.5671 

18.3590 

16.2026 

14.0971 

12.0420 

Center to Center Distance 

(in. ) 

35.0762 

35.6091 

36.0765 

36.4822 

36.8296 

37.1223 

Back Focal Distance 

(in.) 

380 

380 

380 

380 

380 

380 

Graz. Angle at Center of each El. 

(mrad) 

15.0000 

13.5178 

12.0689 

10.6531 

9.2702 

7.9197 

Half Width of Entrance Annulus 

(in.) 

0.24 

0.22 

0.20 

0.175 

0.15 

0.125 

Image Curvature 

(in." 1 ) 

-0.262 

-0.328 

-0.417 

-0.538 

-0.730 

-0.018 

Focal Length 

(in. ) 

398.957 

398.957 

398.957 

398.957 

398.957 

398.957 

Collecting Area 

(in. 2 ) 

71.98 

59.49 

48.30 

37.32 

27.85 

19.83 

Optimized Dist. of Im. PI. 








r Gaussian Focal P. 

for .5 arcsec • 

(10“ 3 in.) 

-8.0 

-8.9 

-9.9 

-11.3 

-12.9 

-15.1 

r 4nl j Max. Half Field 
in Flat Field Ang)e 

(arcmin) 

2.5 

2.4 

2.25 

2.10 

2.00 

1.85 



Image field curvature Is known to be the predominant aberration 
In grazing Incidence telescopes. To illustrate the dependence of the 
curvature on the diameter of the system, the meridional sections of the 
Image surfaces of the Individual systems are shown In Fig. 1.3 In actual 
size. 



Fig. 1.3: Meridional Sections of the Image Surfaces of all 

Individual Subsystems (Actual Size) 
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The field extends over about + 20 arcmin. The resultant optimum fecal 
surface of the compound system Is shown In Fig. 1.4. It Is also drawn 
In actual size and Includes the rms-spot sizes at different field positions 
up to 20 arcmin. Fig. 1.4 was generated by letting each subsystem con- 
tribute according to Its geometric collecting area. The result Is expected 
to change, however, when the effective collecting areas which depend strong 
ly on the grazing angle and the surface material are taken into account. 



Fig. 1.4: Meridional Section of the Optimum Focal Surface Formed by 

all Six Subsystems 
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2. REFLECTIVITY AND EFFECTIVE AREA 


The reflectivities and effective areas {geometric area x reflectivity) 
of the six telescope subsystems were determined for two coating materials 
and four wavelengths. The reflectivities, R, were calculated according to*. 


with 


R * Rj (l+R 0 )/2, 

R = (4a 2 (s-a) 2 + V)/ (4a 2 (s+a) 2 + V), 
1 

R o a (4a 2 (t-a) 2 + V)/ (4a 2 (t+a) 2 + V), 
a 2 = (s 2 -U+ (s 2 -U) 2 + V)/2, 
s = sin y 


t - cos y/tan y, y being the grazing angle. 


The following U - and V - values were provided by L. VanSpeybneck, 


SAO: 

8 

U 

Au 

V 

Ni 

U 

V 

2 

1.54687E-4 

6.21530E-10 

8.41236E-5 

8. 12310E-12 

10.44 

2.57412E-3 

1. 13329E-6 

1.73944E-3 

7.7305E-7 

47.68 

2.39684E-2 

4.38718E-4 

2.95621E-2 

2.27016E-4 

103.32 

1.3276E-1 

1.09378E-3 

7. 16910E-2 

1.08618E-2 


The computed reflectivities and effective areas are summarized on the next 
two pages. 


*B.L. Henke, Phys, Review A, 6, .94(1972) 
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0.079709332 

5.737431332 

2 


0. 01 35 177b 8 

0. 174547064 

1 O . 38330434 

3 


0.012063919 

0. 413903125 

19. 99176245 

4 


9. 010653141 

0 . 605 3 6 3 7 0 2 

22.59217336 

Cj 


9. 27015E-03 

0.705747147 

19. 65505304 

6 


?. 91969E-03 

0.771243543 

15.29375947 




TOTAL EFF AREA : 

93. 65403943 

WAVELENGTH 

< A > : 

19.44 



SYSTEM 


GRAZ ANG 

REFLECTIVITY 

EFF AREA 



(RAD) 


( SQ IN) 

1 


0.01 5 

0.793923223 

57.14659356 

j 


0.013517763 

0.313529435 

48 . 89636606 

3 


0, 012063919 

0. 332741237 

40.22140172 

4 


0. 010653141 

0, 351597144 

:» 1 ■ i o 1 o U 5 4 

5 


9.27015E-03 

0. 370130579 

4 , 3 3 1366 3 



7 , 9 1 9 6 9 E — 0 3 

0. 333370439 

17. 616385 8 




TOTAL EFF AREA • 

9 l 9 . 3959892 

WAVELENGTH 

( A ) : 

47. 68 

■ 


SYSTEM 


u R ft Z A N G 

REFLECTIVITY 

EFF AREA 



RAD) 


• SQ IN) 

1 


0 .015 

0 . 3 865496 3 6 

63.31334640 

2 


0.013517763 

0.397215487 

53 . 37534932 

-J 


0 . 0 1 20639 1 9 

0.907743657 

43.34426014 

4 


0.010653141 

O . 913146716 

34.2652 3 5 4 3 

Cj 


9. 27015E-03 

O. 923406950 

25 . 356 1 3356 

hi 


7.91 969E-03 
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ORIGINAL PAGE IS 
< POOR QUALITY 
3. FOCAL PLANE OPTIMIZATION 

The analyses described in this chapter were performed by tracing rays 
through all six subsystems simultaneously whereby the number of rays for 
each subsystem was proportional to its effective area. 

Since each individual subsystem of the telescope assembly has a 
different focal surface with the field radii increasing with increasing 
aperture diameters, the field radius of the nested system will fall sc > 
where between the radii of the outer and inner systems, defending on the 
contributions of each of the six subsystems to the total effective area. 
Because the subsystem contributions to the collecting area are functions 
of their reflectivities, and thus functions of the wavelength, the field 
radius of the optimum system f . cal surface is expected to be wavelength 
dependent. The best-focus performances and field radii of the nested 
system for two surface materials (Au and Ni) and four different wave- 
lengths are represented by the curves numbered 1 in Figures 3.2 to 3.8. 

To find the optimum position of an assumed flat detector surface, two 
analytical approaches were pursued. 

a) The Common Focal Point Configuration 

The system is aligned so that all six subsystem focal points co- 
incide (See Figure 3.1a). The detector surface is moved toward the 
telescope until the on-axis rms-spot diameter reaches C.5 arcsec. 

The flat-field performances for this focal plane configuration are rep- 
resented by the number 2 curves in Fig. 3.2 to 3.8. The respective de- 
tector shift, A b, from the conmon focal point is also given on the dia- 
gram. 
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b) The Staggered Focal Point Configuration 

The image plane positions for an on-axis rms-spot size of 0.5 
arcsec are determined for each subsystem separately. Matching the six 
image planes then resulted in a staggered focal point array (see Fig. 

3.1 b). The individual image plane shifts from the respective focal 
points are then constant and wavelength independent. These respective 
shifts are from the outer to the inner subsystem: 

Ab = -0.0080 in., 

1 

Ab = -0.0089 in., 

2 

Ab = -0.0099 in., 

3 

Ab = -0.0113 in., 

4 

Ab = -0.0129 in., 

s 

Ab - -0.0151 in.. 

6 

The number 3 curves in Figures 3.2 to 3.8 show the Hat-field performance 
of the staggered focal point configuration for the two metals and the four 
wavelengths. 

Even though the difference between the common focal point configura- 
tion and the staggered focal point configuration seems to be insignificant, 
the wavelength independence of the second configuration must be regarded a 
definite advantage since it avoids the focus change that is required for 
different wavelengths when using the first configuration. 
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A u , 2 i 


1 


HALF FIELD ANGLE (arc mm) 


Fig. 3.2: Off-Axis Performance: (1) Best Focus (Field Radius = -1.65 in.), 

(2) Flat Field, Comnon Foci (Focal Shift = -0.012 in.), (3) Flat 
Field, Staggered Foci. 


A 11,10.44 A 
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HALF FIELD ANGLE (arc min) 


Fig. 3.3: Off-Axis Performance: (1) Best Focus (Field Radius = -2.41 in.), 

(2) Flat Field, Common Foci (Focal Shift = -0.010 in.), (3) Flat 
Field, Staggered Foci. 




RMS-SPOT St2E(«rCMC) 



HALF FIELD ANGLE (arc min) 


Fig. 3.4: Off-Axis Performance: (1) Best Focus (Field Radius = -2.44 in.}, 

(2) Flat Field, Common Foci (Focal Shift = -0.010 in.), (3) Flat 
Field, Staggered Foci. 


o 

e 



HALF FIELD ANGLE Care min) 

Fig. 3.5: Off-Axis Performance: (1) Best Focus (Field Radius = -2.44 in.), 

(2) Flat Field, Common Foci (Focal Shift = -0.010 in.), (3) flat 
Field, Staggered Foci. 
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RMS-SPOT SIZE (arc sec) o> RMS- SPOT SIZE (arc sec) 



HALF FIELD ANGLf (arc min) 


Off-Axis Performance: (1) Best Focus (Field Radius = -1.35 in.), 
(2) Flat Field, Coronon Foci (Focal Shift = -0.013 in.), (3) Flat 
Field, Staggered Foci. 



HALF FIELD ANGLE (arc min) 

Fig. 3.7: Off-Axis Performance: (1) Best Focus (Field Radius = -2.38 in.), 

(2) Flat F'teld, Coircnon Foci (Focal Shift = -0.010 in.), (3) Flat 
Field, Staggered Foci. 
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RMS- SPOT SIZE Carcaac) 



Fig. 3.8: Off-Axis Performance: (1) Best Focus (Field Radius = -2.43 in.), 

(2) Flat Field, Common Foci (Focal Shift = -0.010 in.), (3) Flat 
Field, Staggered Foci. 
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4. EFFECTS OF MISALIGNMENTS ON THE PERFORMANCE 


The effects of misalignments on the telescope performance 
were Investigated In two ways. First the misalignment sensivi- 
tles, l.e., the effects of linear and angular alignment errors 
between primary and secondary on the focal plane performance 
of the Individual subsystems were established. Then the per- 
formance of the entire nested array was determined while each of 
the twelve mirror elements was randomly misaligned within certain 
preset limits. 

4.1 ALIGNMENT SENSITIVITIES 

The following table summarizes the increases of the spot 
size due to various isolated misalignments between the primary 
and secondary mirrors of all subsystems. 

TABLE 4.1: Misalignment Sensitivities 





DE SPACE 

DECENTER 

TILT 

Ui 

2 

UJ 

cm 

o 

f SYSTEM 

1 

0.015 fir ad/ pm 

0.1 firad/ pm 

1 . 9 firad/firad 

SYSTEM 2 

0.013 firad / u m 

0.1 /irad/pm 

1.9 firad/firad 

z 

w 

1x1 4 

SYSTEM 

3 

0.012 /trad/pm 

0.1 fir ad/pm 

1 . 9 firad/firad 

M 
t— t 
U~> 

SYSTEM 4 

0.010 firad/\m 

0.1 /trad/pm 

1.9 firad/firad 

►— 

o 

Cl. 

00 

SYSTEM 

5 

0.09 /irad/pm 

0.1 /irad/pm 

1.9 firad/firad 

00 

i 

.SYSTEM 

6 

0.08 ftrad/pm 

0.1 firad/ pm 

1 . 9 firad/firad 


All spot sizes were established in a fixed gaussian focal plane. 
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4.2 RANDOM ALIGNMENT ERRORS 

The purpose of this experiment Is to predict the most probable 
performance of a real system where each component can only be 
aligned to within certain limits with regard to the perfect de- 
sign. In the computer simulation each of the twelve mirror elements 
was therefore allowed to be out of alignment with respect to five 
degrees of freedom and within given limits. The five degrees of 
freedom were shared by two tilts about the center of each element 
in the x, z - and y, z planes, and linear shifts along the three 
axes. The random misalignment values were generated by multiplying 
a given limit value for every degree of freedom with a computer 
generated random number between -1 and +1. One hundred computer 
runs per set of limit values, each determining the rms spot size 
of the entire telescope system were then made. The performance 
predictions for various misalignments are shown in Figs. 4.1 to 
4.7. While some runs were made where the computer automatically 
searched for t!ie plane of best focus, other analyses were done in 
a fixed gaussian focal plane. The first four graphs show the per- 
formance for angular misalignments of 5 yrad (1 arcsec) and 
linear misal iqnments (decenter and despace) of 0.001 in. In the 
following two runs the tilt and decenter limit values were varied 
successively to appreciate there individual contributions. The 
final plot shows a set of misalignment limits that yield a spot 
size in the order of 1 arcsec. 
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Fig. 4.1: Probable Per- 
formance for the 
Following Limit Values: 
Tilt * 5 yrad 
Decenter 3 0.001 In. 
Despace = 0.001 in. 

(On Axis, Best Focus) 
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Fig. 4.2: Probable Per- 

formance for the 
Following Limit Values: 
Tilt = 5 yrad 
Decenter = 0.001 in. 
Despace = 0.001 in. 

(On Axis, Fixed Focus) 
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Fig. 4.3: Probable Per- 
formance for the 
Following Limit Values: 
Tilt * 5 yrad 
Decenter * 0.001 in. 
Despace = 0.001 in. 

(2 arcmin Off Axis, 

Best Focus) 


PERCENT 
25 

2B 
15 
IB 
5 
■ 



ut 


t ui ‘■I 

- n n' n n‘ t 

RN5 SPOT DIR (RFC StC) 


UI 


Ui 

UI Ul 


Fig. 4.4: Probable Per- 
formance for the 
Following Limit Values: 
Tilt = 5 yrad 
Decenter = 0.001 in. 
Despace = 0.001 in. 

(2 arcmin Off Axis, 
Fixed Focus) 
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Fig. 4.5: Probable Per- 

formance for the 
Following Limit Values: 
Tilt * 5 yrad 
Decenter = 0.005 in. 
Despace =0.001 in. 

(On Axis, Fixed Focus) 
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4.6: Probable Per- 
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28 

Following Limit Values: 
Tilt = 10 yrad 

IE 

Decenter =0.001 in. 
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Despace = 0.001 in. 
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(On Axis, Fixed Focus) 
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PERCENT 


HERN* 1.23MB 
STD. OEV.® 0.2252 


Fig. 4.7: Probable Per- 

formance for the 
Following Limit Values: 
Tilt = 2 yrad 
Decenter - 0.0002 in. 
Despace = 0.0002 in. 

(On Axis, Fixed Focus) 
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5. 


COMMENTS ON TESTING THE X-RAY TELESCOPE 


ORIGINAL PAGE IS 
OF POOR QUALITY 


To perform a meaningful test of the x-ray telescope, it is 
desired to have a point source such that the diameter of the on- 
axis image formed by the telescope does not exceed 1 arcsec. The 
axial spot diameters as a function of the inverse object distance 
for two telescope sizes (1.2m and 1.5m) with equal focal lengths 
of 398.96 in. are plotted in Fig. 5.1. 



l/ OBJECT DISTANCE (in.* 1 ) 


Fig. 5.1: Axial spot size as a function of the object distance 
for two different telescope sizes. 
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IMAGE SHIFT (in. 


The associate focal plane shifts, which are the same for both 
telescope sizes because of the equal focal lengths, are plotted 
in Fig. 5.2 as a function of the inverse object distance. 



1 / OBJECT DISTANCE r in.) 

Fig. 5.2: Fc Shi ft as a Function of the Inverse Object- 

DiSLu.ice for Telescopes With a Focal Length of 
Approximately 400 in. 

There are two principle methods to obtain the required spot 
size. The first and most straightforward method is to select an 
object distance of at least 2000 or 2500 ft., depending on the 
telescope size. The second method would be to use a shorter 
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1QOO ft 


object distance, for Instance, the current vacuun tunnel length of 
1000 ft., and optically project the source to the desired distance. 
In the following, methods for an optical extension of the object 
distance will be discussed. 

a) Use of a single hyperbolic element at a distance of 
1000 ft. and cover the six subsystems of the telescope 
assembly subsequently by changing the source distance 
to the element. 

Comnents: The approach has the same inherent problem 
that the telescope has to begin with, i.e., 
the single element can only be optimized 
for one source distance, and thus for only 
on subsystems of the telescope. Ray trace 
results show that when this distance changes 
in order to cover other subsystems, the spot 
size quickly assumes a diameter of many 
arcsec. 

b) Use a nested system of six single-element hyperbolas to 
simultaneously cover all six subsystems of the x-ray 
telescope assembly. 

Comments: The required diameter of the largest element 

of the test optics is 

D = Do (l - , 
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where Do Is the diameter of the x-ray 
telescope and s Is the desired object dis- 
tance In ft., assuming the test optics Is 
placed at a distance of 1000 ft. Although 
this concept is principally workable, there 
is a problem concerning the required source 
diameter. For a back focal distance, b, 
measured from the center of the test array, 
the source size. Ax, is required to be 
smaller than 5.10~^»b in. in order to obtain 
a spot size of not larger than 1 arcsec. For 
example, a back focal distance of 800 in. 
still requires a source diameter of smaller 
than 100 pm. 
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6. APPENDIX 

KOGRAZ: A RAY TRACE PROGRAM FOR GRAZING INCIDENCE TELESCOPES 

6.1 INTRODUCTION 

KOGRAZ is a computer program that- analyzes grazing incidence 
telescopes based on exact ray tracing. The optical system con- 
sists of two reflective surfaces the shape of which can be any 
conic sections of revolution. The program is written in BASIC 
and performs the following tasks. 

1. Computes x, y-coordinates of incident rays on any 
surface, 

2. Determines centroid coordinates of the image spot 
generated from a point source, 

3. Calculates the rms-image diameter of a point source, 

4. Determines the Radial Energy Distribution (RED) within 
an image spot, 

5. Determines the image-field curvature, 

6. Determines the surface of best focus, 

7. Plots spot diagrams. 

6.2 THE SURFACE EQUATION 

The surface of both elements of two-mirror grazing incidence 
telescopes are conic sections of revolution (see Fig. 6.1). 
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CONIC SECTION 



Fig. 6.1: Meridional Section of Grazing Incidence Surface 

The surface equation is given by 

P 2 - p 2 = 2kz - (1 + S)z 2 , p 2 = x 2 + y 2 

p is the central radius of the element, k is the subnormal at the 
center and 6 is the deformation constant. The type of conic section 
is determined by the value of 5 according to the following list: 


5 > 0 

Oblate Ellipsoid, 

6 = 0 

Sphere, 

0 > 6 > -1 

Prolate Ellipsoid 

6 = -1 

Paraboloid, 

5 <-l 

Hyperboloid. 


Since all elements working in grazing incidence have a defor- 
mation constant very close to -1, the paraboloid, the prolate 
ellipsoid and the hyperboloid are the only types of conic sections 
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of revolution to be considered. 

It may be Interesting to note that the subnormal, k for P o *0, 
i.e. , in the vertex equation, is identical with the radius of 
curvature at the vertex. 

6.3 THE TWO-MIRROR GRAZING INCIDENCE TELESCOPE 

The two-mirror grazing incidence telescope as shown in Fig. 

6.2 is completely defined by the following four parameters: 

Grazing Angle: y (G) 

Center Radius of First Mirror: (R (1) ) 

Center to Center Separation: d (D (3) } 

Half Widths of Entrance Annulus: Ap (T) 

PARABOLOID 



Fig. 6.2: Schematic of Two-Mirror Grazing Incidence Telescope 
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The bracketed expressions are the corresponding Input parameters 
used in the program. All other quantities can be expressed by 
using only the first three input parameters. The following is 
a simnary of the most useful system parameters and the relations 
among them. 

Center Radius of Second Mirror: 

Back Focal Distance (measured 
from center of second mirror): 

Center Subnormal of First 
Surface: 

Center Subnormal of Second 
Surface: 


Deformation Constant of First 
Surface: 

Deformation Constant of 
Second Surface: 


p o2* p o1' d * tan2Y 
b = P o2 /tan 4y 

k l * - p ol tan Y 
k 2 * - p q2 tan 3 y 

= 'I 

, = _ | sinfr ] 

2 [ sin 4 y - sin 2y | 


System Focal Length: f = 2p o -j/sin8y 

System Focal Ratio: F = 1/2 tan 4y 

Total Length of both . = 2 a 

Elements: 
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Entrance Aperture: The width of the entrance annulus is calculated 

by inserting ^uii\AL PAGE ifa 
OF POOR QUALITY 

Z = -d/2 and Z = + d/2 
into the equation of the first surface, 

P 2 - P 2 , ■ k, z. 

a) p = d q1 + Ap-j and k 1 = yo ol . 

2 2 2 

It then follows: p ^ + 2 p q1 Ap^ + Ap-j - p - = 2o ol Ap-j = 2yp , d/2, 

or Ap-j = y d/2 . 

b) p = p j - Ap 2 , and - -yp ^ 

It then follows: Ap^ = yd/2. 

The half widths of the entrance annulus therefore is to a good 
approximation, 

Ap = Ap^ = Ap 2 = yd/2. 

The collecting area is given by 

A = tt(p 0 -| + Ap ) 2 - tt(p o1 - Ap ) 2 

= 4 it p i Ap = 2 it p ,yd. 
oi or 
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6.4 INPUT PARAMETERS AND PROCEDURE 


The input consists of parameters describing the optical system 
and the initial parameters defining the set of rays traced through 
it. The optical system is characterized by its surfaces, their 
positions and orientations, and the entrance aperture. 

A sinroary of all pertinent input parameters is listed below. 



Fig. 6.3: Ray Coordinates in Entrance Aperture 

Fig. 6.3 explains the ray distribution in the entrance annulus. 
T is the radial ray separation and A is the angular tangential ray 
separation. A is generally calculated so that the linear tangential 
ray separation equals T. 
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TABLE 6.1: Input Parameters 


Name 


Description 


Line-No. 


J Total number of surfaces, (incl. dummy surfaces) 

G Grazing angle at center of each surface (rad) 

GO Off-axis angle in XZ-plane, (rad) 

HO Off-azis angle in YZ-plane, (rad) 

SO Divider of T 

T Half width of entrance annulus 

A Ray input parameter, (rad), (See Fig. 3) 

MO Inverse object distance (M0 = 0 for telescope) 

PI Starting point for RED in percent 

P 0 Increment of RED in percent 

R(l) Center Radius of first surface 

D(l) Distance from entrance aperture to center of 

first surface 

D{3) Mirror center to center separation 


10 

40 

50 

55 

60 

70 

80 

90 

95 

96 
170 

180 

190 


F 


S(3) 
T(3) 
G(3) 
H ( 3 ) 
C ( 3 ) 
C(5) 


Unit selector for image size 


’ F = 1 {linear units] 

F = Focal length [rad] 450 

F=J0‘Focal length/3600 
[arcsec] 


Tilt of secondary about X-axis, (rad) " 

Tilt of secondary about Y-axis, (rad) 

Decenter of secondary in x-di recti on 460 

Decenter of secondary in Y-di recti on ” -490 

Despace primary - secondary 
Defocus 
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Finally there are the subroutines which are those parts 
of the program that perform the tasks listed In the Introduction. 
After loading the program, all subroutines are bypassed using 
a "GOTO" statement. A temporary elimination of this statement 
..lakes the corresponding subroutine available. The following 
listing summarizes all possible print-out options: 


TABLE 6.2: Print Out Options 


Parameter listing 

Ray coordinates in entrance plane and on J-th 
surface 

Best-Focus computation 

Centroid coordinates 
Field-curvature, RMS-Spotsize 

Radial Energy Distribution 

Percent of rays through specified exit aperture 


300 GOTO 450 

1380 

GOTO 

1440 

1580 

GOTO 

1740 

1760 

GOTO 

1910 

1920 

GOTO 

2240 

2270 

GOTO 

2300 
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ORIGINAL PAGE IS 
OF POOR QUALITY 
6.5 PROGRAM LISTING 


1 REM **»***»»*»*#»****«#»KOGRAZ*****»************** 

10 J = 5 

20 J0=ATN( 1)/45 
30 F=J0*399/3600 
40 G=0 .0 15 
50 G0=0 
55 H0=0 
60 SQ = 2 
70 T=0.245 
60 A= 16 * JO 
90 MOsO 

95 P 1 =60 

96 P0=10 

ICO DIM U(500),V(500),X(500),Y(500) 

110 DIM A(500) ,B(500) ,R(5C0) ,Z(500) 

120 FOR N= 1 TO J 

130 R ( N ) =0 , L( N ) =0 , D( N ) = 0 , G ( N ) = 0 , K ( N ) =0 ,S ( N ) =0 , T ( U) =0 ,C ( N ) =0 

140 K(N)=1 

150 E ( N ) = - 1 

160 NEXT N 

170 R C 1 ) = 2 3 • 6 8 

160 D( 1 )= 16.5 

190 D(3)=35. 076156 

200 K ( 1 ) =- 1 , K ( 3 ) =- 1 

210 A 1 =2*G 


220 A2=4*G 

2 30 R ( 3 ) = R ( 1 ) -D ( 3 ) *TA N ( A 1 ) 

240 D(5)=H(3)/TAN(A2) 

2 50 £(3)=-(SIN(A2-A1 )/ (SIN ( A2) -SIM ( A 1 ) ))*2 
260 L( 1)=-R( 1)*TAN(A1/2) 

270 L(3)=-R(3)*TAN(A1/2+A2/2) 

275 F=2»R: 1)/SIN(6*G) 

260 IF G0+H0=0 THEN 290 

264 PRINT 

265 PRINT "FIELD ANG : GO="GO, "H0="H0 

266 PRINT 


290 REM ********** DESIGN PARAMETERS********** 
300 GOTO 450 
330 PRINT 

340 PRINT "GRAZ ANG :";G 
350 PRINT 

360 PRINT "RAD OF EL :";R(1),R(7) 

390 PRINT 

400 PRINT "SEPARATION: " ;D( 1 ) ,D( 3) , D(5) 

410 PRINT 


420 PRINT "K-CONSTANT: " ; L( 1 ) ,L( 3) 

4 3C PRINT 

^40 PRINT "D-CONSTANT : " ; E ( 1 ) , E ( 7 ) 

441 PRINT 

445 PRINT "FOCAL LENGTH: "F, » T:"T 

4 50 F = F 
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500 «3 = 0,MsO,iMsO 

510 FOR H 0 = R ( 1 ) — T TO R(1)+T STEP T 
520 FOR B = A TO 360*J0 STEP A 
530 XO=RO*COS(B) 

540 Y0=R0»SIN(B) 

550 X=X0 
560 Y=Y0 
570 Z = 0 

580 K5=TAN(G0) -MO*X 
590 L5=TAN(H0)-M0*Y 
600 C0=1/SQR(HK5~2+L5~2) 

610 A0=K5*C0 
620 B0=L5*C0 
630 FOR 1=1 TO J 
640 D1=C(I) 

650 D2=H(I) 

660 D3=D(I)+C(I) 

670 W1sS(I) 

680 rt2=T ( I ) 

690 n = L ( I ) 

700 D=E ( I ) 

710 K=K(I) 

720 IF w 1 *W 1+W2*W2+W3*w3=0 THEN 630 
730 A1=COS(W2)*COS(W3) 

7 40 B 1 =CQS ( W 1 ) * S I N ( W 3 , +S I N ( W 1 ) * S I N ( W 2 ) *COS ( w 3 ) 
750 C 1 sSIN ( W 1 ) «SIN ( W3 ) -COS ( W 1 ) *51 N ( w2 ) *COS ( w 3 ) 
760 A2=-CCS(W2)*SIN(W3) 

7 70 B2=COS( W 1 ) *COS ( Wj ) -SI N ( In 1 ) *SI N ( Vv 2 ) *SI N ( W 3 ) 

7 cG C 2 = S I N ( w 1 )»C0S(«3)+C0S(« 1 ) *SI N ( '*2 ) *SIN f W 3 ) 
790 A3=SIN(*2) 

800 B3=-SIN(rt1}»COS(W2) 

810 C3=COS(W1)*COS(W2) 

620 GOTO 850 
c 30 A 1 = 1 , 5 2 = 1 , C 3 = 1 

8 4 0 A 2 = 0, A3 = 0,B1=0,B3 = 0,C 1=0,02 = 0 
650 X1=A1*(X-D1)+B1*(Y-D2)+C1*(Z-03) 

660 Y1=A2*(X-Dl)+b2*(Y-D2)+C2»(Z-D3) 

670 Z1rA3*(X-D1 )+B3*(Y-D2)+C3 # (Z-D3) 

000 A4=AG*A1+BO*B1+CO*C1 
090 54=A0*A2+B0*52+C0*C2 
900 C4=A0»A3+B0*d3+CJ*C3 
910 IF R = 0 THEN 1020 
930 51=D+1/C4~2 

9 40 32 = R-A4*(X1-A4*Z!/C4)/C4-rj4’ , CY1-B4*Z1/C4)/C4 
950 S3=(X1-A4*Z1/C4)*2+(Y'-B43Z1/C4)'2-R(I)"2 

10 00 Z=S3/(S2*( USQR( 1-S1*S3/o2“2) ) ) 
loio core 1030 

1020 z=o 

1030 X=A4*(Z-Z1 )/C4+X1 
1040 Y=B4*(Z-Z1)/C4+Y1 
1050 IF 5=0 THEN 1110 
1060 P = fi-(UD)*Z 
1070 Z 2 = X / P 
1090 Z 2 = x / P 

1 ICO GOTO 1 125 
1110 Z2=0 

1 120 Z3 = 0 

1125 Z4 = GQH(Z2*Z2 + Z2*Z3+1 ) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

1130 A5=Z2/Z4 

1 140 B5 = Z3/Z4 

1 150 C 5 = - 1 / Z 'J 

1160 V=A4*A5+B4*B5+C4*C5 

1200 U=V*( 1-KJ/X 

1220 A0=A4/K-U*A5 

1230 B0 = 64/K.-';»B f ; 

! 240 C0=C4/K-u*C5 

1250 IF IOI ThEN 1270 

1260 M = M+ 1 

1270 NEXT I 

1280 N=N+ 1 

1 290 L) ( N ) = XO 

1300 V(N)=YO 

1310 X(N)=X 

1320 Y(N)=Y 

1330 A(N)=A0/C0 

1340 B (N ) =B0/C0 

1350 NEXT B 

1350 NEXT RO 


1370 REM ***»**«*»*raY COORDINATE ;»»»******* 

1380 GOTO 1440 

1390 PRINT " XO YO X Y" 

1400 PRINT 


1410 FOR J=1 TO N 

1420 PRINT U ( J ) ,V(J) ,X(J) ,Y(J) 

1430 NEXT J 
1440 PRINT 

1450 X9=0, Y?=0, A9=0,B9=0 
1460 FOR J=1 TO N 
1470 X9=X9+X(J) 

1460 Y9=Y9+Y(J) 

1490 A9=A9+A(J) 

1500 B9=B9+BfJ) 

1510 NEX1 J 
1520 X 9 = X 9 / N 
1530 Y9=Y9/N 
1540 A9=A9/N 
1550 B9=E9/N 
1560 Z9=0 

1570 REM »*«****»**3EST FOCUS***** ***** 

1580 GOTO 1740 
1590 Z7=0,Z6=0 
1600 FOR J=1 TO N 

1610 Z7=Z7+(X(J)-X9)*(A(J)-A9)+(Y(J)-Y9)*(B(J)-B9) 

1620 Z8 = Z8 + (A(J,'-A9)~2+(5(J)-B9)~2 

1630 NEXT J 

1640 Z9=-Z7/Zo 

1650 X9=0, Y9=0 

1660 FOR J=1 TO N 

1670 X(J)=X( J)+A(J)*Z9 

1 6 G 0 Y( J ) =Y( J ) +B( J ) *Z9 

1690 X9=X9+X( J) 

1700 Y9=Y9+Y(J) 

1710 NEXTJ 
1720 X 9 = X 9 / N 
1730 Y Q = Y 9 / N 
1740 REM 



1750 HEM *»»****»* *C ENT HO ID FIELL-CURV RMS-SPOT SIZE 
1760 GOTO 1910 
1770 R6 = 0 
1760 FOR J=1 TO N 

1790 R6=m 6+(X(J)-X9)~2+(Y(J)-Y9)~2 
1600 NEXT J 
1010 H9=SOh(R6/N) 

1620 PRINT "CENTROID X=";X9 
1630 PRINT "CENTROID Y = " ; Y5 
1640 PRINT " DELTA Z=";Z9 
1650 PRINT 
1 6b0 IF GG~2+h0~2=0 ThEN I8b0 
1670 PRINT "FIELD-CURV=";2*Z9/(X9*2+Y9“2) 

1 660 PRINT 
1690 R9=R9/F 

1900 PRINT "RMS-DIA :";2*R9 
1910 PRINT 
1920 COTO 224G 
1925 FOR G=1 TO N 

1930 Z(J)=SCT ((X( J)-X9)*2+(Y(J)-Y9)“2) 

1 940 NEXT J 
1950 FOR K= 1 TO N 
I960 D=Z( 1 ) 

1970 H1 = 0 

1930 FOn J=2 TO N-K+1 
1990 IF Z ( J ) >D TnEN 2030 
2000 Z(J-1 )=D 
2C1C D=Z(J) 

2C2C GoTO 2060 
2030 E=Z ( J ) 

2G4C Z(J)=Z(J-1) 

2C5C Z ( J - 1 ) = E 
2CoO NEXT J 
2070 R(K)=D 
2060 NEXT K 

2090 REM **********RADIAL ENERGY DISTRIBUTION********** 

2110 PRINT "RADIAL ENERGY DISTRIBUTION:" 

2120 PRINT 

2130 PRINT " % DIA" 

2 140 PRINT 
2150 11=0 

2160 FOR 10 = F 1 TO 100 STEP PO 

2170 FCR 1=11+1 TO N 

2 IcG IF I* 100/N<10 ThEN 2220 

2 190 PRINT 10, 2*R ( i)/F 

2200 11=1 

22 1 C GOTO 2230 

2220 NEXT I 

2230 NEXT 10 

22^0 PRINT 

22 50 PRINT "NO. uP RAYS: -, ;N 

2260 REM *»****»*«« PE, nCENT OF RAYS THROUGH EXIT APERTURE********** 
2270 GuTO 2300 
2 2o0 PRINT 

2290 PRINT N* 100/M ; " 5 ThROUCh EXIT APERTURE" 

2300 PRINT 
9999 EuD 
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